Self-consistent relativistic random phase approximation 
with vacuum polarization 
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We present a theoretical formulation for the description of nuclear excitations within the frame- 
work of relativistic random-phase approximation whereby the vacuum polarization arising from 
nucleon-antinucleon fields is duly accounted for. The vacuum contribution to Lagrangian is explic- 
itly described as extra new terms of interacting mesons by means of the derivative expansion of 
the effective action. It is shown that the self-consistent calculation yields zero eigenvalue for the 
spurious isoscalar-dipole state and also conserves the vector-current density. 
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INTRODUCTION 

The relativistic field theory based on the quan- 
tum hadrodynamics (QHD)[l| has been very suc- 
cessful in describing the nuclear properties not 
only for ground states but also for excited states. 
Although the response of a system to an external 
field has been already investigated in the eighties 
by using the relativistic random-phase approxima- 
tion (RRPA) in the relativistic mean-field (RMF) 
basis H H H H H H i|, the self-consistent 
methods with nonlinear effective Lagrangian for 
a quantitative description of excited states have 
been developed only during the last few years 
H EH IH El However, in particular, it is 
important to emphasize here that the negative- 
energy RMF states contribute essentially to the 
current conservation of RRPA eigenstates and the 
decoupling of the spurious state. In our recent 
study, it has been shown that the negative-energy 
RRPA eigenstates generated from the RRPA equa- 
tion with the fully consistent basis have a signif- 
icant role for gauge invariance in the electromag- 
netic response|l4). 

Almost all investigations of RRPA in the QHD 
model referred to above, however, neglected actual 
antinuclcon degrees of freedom. The basis set used 
for the RRPA calculation is usually obtained from 
the RMF theory wherein only positive-energy nu- 
cleons are taken into account and the Dirac sea is 
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always regarded as unoccupied (this is sometimes 
called no-sea approximation). This approximation 
is very convenient because we do not have to puzzle 
over a renormalization procedure not only in the 
calculation of basis set under the full one-nucleon- 
loop contribution, which we refer to as the rela- 
tivistic Hartree approximation (RHA), but also in 
the RRPA calculation in which the Feynman part 
is essentially divergent. Thus, the RRPA calcu- 
lations in a finite nuclear system with the inclu- 
sion of vacuum polarization have been performed 
only approximately 

HE 110 Au th ° RHA + 

RRPA calculations performed earlier have been 
based on the local-density approximation, that is, 
the renormalization of the one-nucleon loop in the 
RHA calculation and of the Feynman part in the 
polarization function were done in nuclear mat- 
ter and the results were applied to finite nuclei. 
However, as the local-density approximation of 
the nucleon-antinucleon loop corrections violates 
the self-consistency for finite nuclei, the spurious 
isoscalar-dipole strength associated with the uni- 
form translation of the center-of-mass dose not 
get shifted all the way down to zero excitation 
energy @. 

The main aim of this letter is to demonstrate 
as to how such a deficiency contained in the pre- 
vious RHA + RRPA calculations is removed by 
employing the derivative-expansion method to es- 
timate the vacuum contribution. We verify the 
consistency by calculating explicitly the spurious 
state in the isoscalar-dipole mode and the current 
conservation of transition densities. 

In the following we work out the vacuum- 
polarization effects using the Lagrangian density 
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of the Walecka a — lo model which is given by 



effective Lagrangian density is given by 
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-(dfj.uju - duuo^) 2 + -m^uo^ - SC, 



(1) 



where U{a) = grS^c 3 + \\gz°~ A denotes the self- 
interaction terms of scalar meson, and SC = 

represents the counterterms to regularize the nu- 
cleon self-energy. The renormalization procedure 
in a finite nuclear system requires considerable ef- 
fort even at the mean- field levelUE 113- 

The ex- 
plicit calculation of the vacuum polarization for fi- 
nite system has been performed recently by Haga 
ct. al. 15] wherein the density variation has been 
found to be substantially large. The use of the ef- 
fective action developed in Ref. [l7|. on the other 
hand, provides a direct and simple approach to es- 
timate the vacuum corrections. It is interesting 
to observe that the lowest order of derivative ex- 
pansion for the one-baryon-loop correction agrees 
with the rigorous results [lij. In the derivative ex- 
pansion, the contribution from the Dirac sea to 
the Lagrangian density, which is expressed by the 
trace and the logarithm of the inverse of the Dirac 
Green function, is expanded by derivatives of the 
meson fields as expressed by 



d A p 

——[Trln(-/^ - m N + g a a - g^Lo^) 
( 27T ) 
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= -V F {a) + -Z° F {a){d,aY 

+^Z u F {a){d ll w v -d v u il f + --- . (2) 

Each term of the right-hand side of Eq.(J2} is fi- 
nite, and therefore we can treat them as ordi- 
nary potential terms. Here, we may need to 
use a standard technique of the renormalization 
to obtain explicit forms of Vf(u), Z F (a), and 
Z F (tj). The method of calculatin g t hem has 
been discussed by many authors [Til ll^. Eof and 
it has also been verified that the convergence of 
the expansion is quite rapid within a mean-field 
aDProximation flil lid . Elj . In the present calcula- 
tion, we employ only the first three terms in the 
right-hand side of Eq. J5J . 

We proceed now to describe the calculational 
details of the random-phase approximation based 
on the leading-order terms of the derivative ex- 
pansion. In the a — lo model, we know that the 
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where superscript + in the nucleon-field operators 
means that only the positive-energy states are to 
be treated explicitly. Then, assuming the station- 
ary and spherical system, the ground state expec- 
tation values of the a and loq, which are written as 
4> and Vo, satisfy the following coupled equations: 

+ m 2 )0 = 5ct <^ + ) - U'tf) - V F {4>) 



-Z F '(<j>)(d^) 2 -d^Z^d^) 



+ -z-'(4>)(d,v,)\ 
+ d li (z F -(4>)d»v ), 



(4) 
(5) 



where the contributions of the negative-energy 
states to source terms are contained in Vf and 
Z F . The first step in calculating the RRPA re- 
sponse is the computation of the RHA, which is 
just to solve these equations together with the 
Dirac equation under the potentials of <f> and Vq. 
The potentials achieved in the mean-field calcu- 
lation are then used to generate the lowest-order 
polarization function. We mention here that the 
Vf and Z F terms, which appear in the meson 
propagator parts, correspond to the excitation of 
negative-energy states to all the positive-energy 
states. Thus, we solve the RRPA equation given 

by 

n RPA (T a ,T b ;p,q;E) = Il D (T a ,T b ;p,q;E) (6) 

+Y,9i f dk 1 dk 2 n D (r a ,T i ; P ,k 1 ;E) 

x D l (k 1 ,k 2 ;E)U RPA (r,T b - fc 2 , q; E), 

where the summation is over the meson fields and 
T's are the 4x4 matrices which denote the ver- 
tex couplings. This RRPA equation exactly has 
the same form as that of no-sea approximation. 
However, it should be noted that there is an es- 
sential difference in regard to the employed me- 
son propagator, Di, in which we must include the 
vacuum-polarization corrections to perform the 
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TABLE I: Parameter set used in the present work. 



m CT (MeV) 


m^(MeV) 


g a g 2 (fm 1 ) 


53 


RHA 458.0 


814.0 


7.10 24.09 


-15.99 8.85 



consistent calculation. The a — u> coupling term 
Zp(a)(d fl ui l/ — dyUJ^) 2 in the present Lagrangian 
provides the coupled equations between the a- 
meson propagator and the time component of the 
w-meson propagator 
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MeV, respectively as the function of the excita- 
tion energy. The most remarkable result of our 
RRPA calculations including the vacuum polariza- 
tion is that the spurious state which is the collec- 
tive mode corresponding to the center-of-mass mo- 
tion, appears at zero excitation energy in both nu- 
clei (shown by solid curves) as also obtained in the 
conventional RRPA calculation without vacuum 
polarization (shown by dashed curves) 0. The 
spurious state is clearly separated from the phys- 
ical states, which are shown by using the small 
imaginary part of energy rj = 0.05 MeV. We em- 
phasize that although earlier RRPA calculations 
with vacuum polarization have never succeeded to 
decouple the spurious state 0], we are now able to 
achieve this by handling the Lagrangian J3J cor- 
rectly. 



while the spatial component of w-meson propa- 
gator can be evaluated independently. Here, the 
Fourier transforms of 
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and the effective meson propagators Di as well as 
free meson propagators D® are expressed as the 
matrices of the momentum space. Details of the 
formulation on the present technique along with 
the comparison between the Feynman part in our 
calculation and that in the local-density approxi- 
mation will be elucidated in a forthcoming publi- 
cation. 

The parameter set used in the present work is 
listed in Table I and has been determined to repro- 
duce the total binding energies and charge radii 
of the spherical nuclei in the RHA calculation. 
This parameter set is similar to that introduced in 
Ref. [22| , where the derivative expansion has been 
used. Small difference between the two sets stems 
from the fact that in the treatment of Ref. the 
one-meson loop corrections have been considered 
in the calculation of the ground state. 

In what follows we discuss the results of our 
RRPA calculations. In Figs. 1(a) and 1(b) we dis- 
play the distributions for the Coulomb responses 
of isoscalar-dipole mode in 16 O and 40 Ca at the 
momentum transfers of q — 237 MeV and q = 118 
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FIG. 1: Distributions of isoscalar-dipole strength in 
(a) 16 O and in (b) 40 Ca. 

Another important aspect for the correctness of 
the RPA calculations is to check whether transi- 
tion charge density (I'\\M\(q)\\I) and current den- 
sity (7'||T\L(g)||/), connecting the ground state I 
and the excited states I' for different multipolar- 
ity A and L, satisfy the conservation law 0, 0|- 
Assuming lun denotes the excitation energy of the 
nucleus, the conservation relation is given by, 



uj N (I'\\M x (q)\\I) 
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<nl ^±(I'\\T xx+1 (q)\\I). 



(12) 
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FIG. 2: RHA + RRPA results for the current con- 
servation in the transition densities for the lowest 
positive-energy (12.94MeV) isoscalar-dipole state in 
16 0. The A{) depicted by thick-solid line shows the 
extent of violation of current conservation. For details 
refer to the text. 

Our results for the RRPA transition current have 
been depicted in Fig. 2 where the contributions 
from the left-hand side (Ihs) and right-hand side 
(rhs) of Eq. I|12f> are shown separately by the 
dashed and dash-dotted lines, respectively. Note 
that for the purpose of clarity the Ihs of Eq. ljT2l 
has been plotted with a negative sign. The differ- 
ence of these two contributions denoted by Ap and 
shown by thick-solid line in Fig. 2, represents the 
violation of the current conservation. It is grati- 
fying to see from Fig. 2 that the RRPA transition 
current is sufficiently conserved. 

In summary, we have studied the self-consistent 
RHA + RRPA method including the vacuum- 
polarization contribution given by the derivative- 
expansion method. We stress here that the present 
method is able to treat the change of the negative- 
energy states due to the presence of particles in the 
positive-energy states for the formation of the nu- 
cleus and also to treat consistently excited states 
using the RRPA. In contrast to the previous calcu- 
lations based on the local-density approximation, 
the self-consistency of the model has been fulfilled 
so that we have obtained the desired results for 
decoupling the spurious isoscalar-dipole state, and 
for conserving the current density. 

Also, it would be pertinent to take stock of the 
present scenario as regards to where we stand now. 
Indeed we now have a powerful method of de- 
scribing the ground state of nuclei in the relativis- 
tic mean-field approximation with the inclusion of 
negative-energy states which are influenced by the 
meson mean fields. Further, we also now have a 
method of calculating the excited states in nuclei 
whereby the excitation of nucleons from negative- 
energy states to positive-energy states is duly in- 
cluded. This makes the present treatment of ex- 
cited states fully consistent with the description 



of the ground state. Numerically, however, the a- 
and w-mean fields obtained in the present treat- 
ment are found to be insufficient in strength, en- 
tailing a weaker spin-orbit splitting which is about 
a half as compared to that obtained in the usual 
RMF model |24| . There are several possibilities to 
provide the missing spin-orbit splitting. One inter- 
esting idea is to explore the possibility of surface 
pion condensation suggested recently j^. An- 
other thing could be to consider the tensor cou- 
pling for the u> meson (2^, E3| ■ These extensions 
of the model constitute an effective field theory 
including the vacuum polarization. Further cal- 
culations testing the RRPA with these extended 
RHA models for a wide variety of observables in 
nuclear excitations would be of immense interest. 
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